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ABSTRACT 

A large fraction of the gas in galactic haloes has temperatures between 10"*'^ and lO'' K. 
At these temperatures, cooling is dominated by metal-line emission if the metallicity 
Z > 0.1 Zq and several lines may be detectable with current and upcoming instru- 
ments. We explore this possibility using several large cosmological, hydrodynamical 
simulations from the Overwhelmingly Large Simulations project. We stack surface 
brightness maps centred on galaxies to calculate the expected mean surface brightness 
profiles for different halo masses. We only consider emission from gas with densities 
riH < 0.1 cm^"^. Results are shown for soft X-ray lines at z = 0.125, ultra-violet (UV) 
lines and hydrogen Balmer-a (Ha) at z = 0.25, and rest-frame UV lines at z = 3. 
Assuming a detection limit of 10~^ photon s^^ cm~^ sr~^, proposed X-ray telescopes 
can detect Oviii emission from z — 0.125 out to 80 per cent of the virial radius 
(i?vir) of groups and clusters and out to 0.4i?vir for haloes with masses 10^^~^^ Mq. 
Emission lines from Cvi, Nvii, Ovil, and Ne X can be detected out to smaller radii, 
0.1 — 0.5i?vir- With a detection limit of 10~^° erg s~^ cm~^ arcsec"^, future UV tele- 
scopes can detect Cm emission out to 0.3-0.6i?vir at z = 0.25, depending on halo 
mass. C IV, O vi. Si iii, and Si iv can be seen out to 10 — 20 per cent of the virial radius 
in haloes more massive than 10"'^^ Mq. Optical Ha emission is comparable in strength 
to Cm emission. At z = 3 it may be possible to observe Cm out to 0.2 — 0.3i?vir 
and other rest- frame UV lines out to ~ 0.1i?vir in haloes larger than 10^^ Mq with 
upcoming optical instruments. Metal-line emission is typically biased towards high 
density and metallicity and towards the temperature at which the emissivity curve 
of the corresponding metal line peaks. This bias varies with radius, halo mass, and 
redshift. The bias is similar for the different soft X-ray lines considered, whereas it 
varies strongly between different UV lines. Active galactic nucleus (AGN) feedback 
can change the inner surface brightness profiles significantly, but it generally does not 
change the radius out to which the emission can be observed. Metal-line emission is 
a promising probe of the warm and hot, enriched gas around galaxies and provides a 
unique window into the interactions between galaxies and their gaseous haloes. 

Key words: galaxies: evolution - galaxies: formation - galaxies: haloes - diffuse 
radiation - intergalactic medium ~ cosmology: theory 



1 INTRODUCTION 

Haloes grow by accreting gas from their surroundings, the 
intergalactic medium (IGM), which is the main reservoir of 
baryons. Galaxies grow by accreting gas from their haloes, 
from which they can form stars. Some of the gas is re- 
turned to the circumgalactic medium (CGM) in galactic 
winds driven by supernovae or active galactic nuclei (AGN). 
Metals produced in stars are also blown out of the galaxy 
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and enrich the CGM. Hence, to understand the evolution of 
galaxies, one needs to study the evolution of the halo gas. 

Ultra-violet (UV) and X-ray absorption line studies 
have revealed both the cold, neutral gas and the warm-hot 
intergalactic medium around galaxies. Unfortunately, this 
type of observation can only provide information about the 
gas along the line of sight to the background quasar, so there 
is no information about the transverse extent of the absorb- 
ing gas cloud. Another limitation is that we can only probe 
gas in case there is a bright background object. Therefore, 
although absorption line studies are excellent for studying 
regions of moderate overdensity in the IGM, the statistics 
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are poorer for relatively dense halo gas, particularly around 
massive haloes. 

The gas emissivity scales with the square of the den- 
sity. Emission is thus dominated by high-density regions 
and is thus a better tracer of the CGM than of the gen- 
eral IGM. Emission line studies thus complement absorption 
line studies. They have the added advantage of providing a 
two-dimensional image, in addition to the third dimension 
provided by the redshift of the emission line, allowing us to 
study the three-dimensional spatial distribution. 

H I Lyman-Qf (Lya) emission originates from regions 
with temperature T ~ 10'* K. It therefore provides an 
excellent route to studying cold gas in the CGM. Diffuse 
Lya emission has been detected ( statistically) around high- 
redshift galaxies at g ^ 3 (e.g. ISteidel et all I2OO0I . I2OIII : 
iMatsuda et al.1 12004 iBower et al.l 120041 ). However, depend- 
ing on halo mass, a large fraction of the halo gas is ex- 
pected to heat to temperatures above 10*'* K, either through 
photo-ionization, accreti on shocks, or through shocks c aused 
by galactic winds (e.g. Ivan de Voort fc Schavd 12012! ) . UV 
metal-line emission enables us to probe gas with T — 



10* 



K, which tends to be more diffuse than colder Lya 



emitting gas and thus a better probe of the average gas 
properties around galaxies (e.g. bertone fc Schave 2013 ). 
Soft X-ray metal-line emission traces even hotter gas with 
T = 10^~^ K, which is the dominant temperature of the 
halo gas around high-mass galaxies and in galaxy groups. 
At T ~ lO'* ''-'' K the emission is do minated by metal lines 
for ^ > 0.1 Zq (e.g lWiersma et al.i r2009a) rather than by 
hydrogen lines or continuum emission as is the case for lower 
and higher temperatures, respectively. 

Dilute halo gas that has been shock-heated to the 
virial temperature is routinely detected in X-ray obser- 
vations of the centres of clusters and groups of galaxies, 
and it may even have been s een around individual galax- 
(e.g. ICrain et all l2010al lbl: I Anderson fc BregmanI I20T1I : 



iBertone fc Schavell2012l : iRoncarelli eraLll2012l ) . These stud- 
ies are focussed mostly on quantifying the emission expected 
from the warm-hot intergalactic medium in a large section 
of the Universe or in mock datacubes. 

It is possible to increase the signal-to-noise by stack- 
ing many observations, centred on galaxies. In this way we 
cannot study the halo gas around a single object, but we 
can characterize the general properties of gas around a cer- 
tain type of galaxy. In this work we will use cosmological, 
hydrodynamical simulations from the overwhelmingly large 
simulations (OWLS) project (Schave et al. 2010.) to quantify 
the expected surface brightness (SB) for the brightest metal 
lines for a range of halo masses, in the (rest-frame) UV at 
high {z = 3) and low {z = 0.25) redshift and in soft X-ray at 
low (2 — 0.125) redshift. Our study complements previous 
work on metal-line emission from the IGM by focussing on 
the CGM and by predicting mean surface brightness profiles 
for different halo masses. 

This paper is organized as follows. In Section [2] we de- 
scribe the simulations we used, how we identify haloes, and 
how we calculate the emission signal. The detection limits of 
soft X-ray lines at z — 0.125 and UV lines at 2 = 0.25 and 
2 = 3 with current and future instruments are described 
in Section |3l The obtained surface brightness profiles are 
shown and discussed in Section [l] Section [5] discusses how 
the gas temperature, density, and metallicity of gas traced 
by metal-line emission are biased. The effect of AGN feed- 
back on our results is discussed in Section [6] We conclude in 
Section[7]and resolution tests are discussed in the appendix. 

In all our calculations we assume the same cosmological 
parameter values as were used during our simulations (see 
Section O = 1 - Qa = 0.238 and h = 0.73). 



2 METHOD 



iLi fc Wang|l2012l ). Most of the detectio ns are made in X-ray 2.1 Simulations 



continuum emission and iron lines, see iBohringer fc Werneii 
l|2010l) for a recent review, but many other X-ray lines 
have been identified. In the last year, a lot of progress has 
been made on obs erving X-ray e mission cl ose to the virial 
radii of clusters (ISimionescu et al. 2 01 li: Akamatsu et al.l 
2OIII: iMiller et all I2OI2I : lUrban elTall l201ll) . but see also 



Eckert et al.l ( 201 ll ). Additionally, there are claims of de- 



tections of metal-lin e emission from the warm-hot inter- 
galactic medium (e. g. Kaastra et al. |200^: iTakei et al.ll2007l : 



IWerner et all I2OO8I : iGaleazzi et all |2012| ). whereas other 
observations set upper limits at lower valu es than these 
claimed detections (e.g. lMitsuishi et al ] |2012l . and references 
therein) . 

Observing metal-line emission from diffuse halo gas 
would yield a lot of information about the physical state 
of the gas, the distribution of metals, and thus the cycle of 
gas between haloes and galaxies. In general, however, emis- 
sion from gas outside of galaxies is faint and thus difficult 
to detect. Many missions have been proposed to study the 
diffuse halo gas in emission. The next generation of spectro- 
graphs should be able to detect certain metal lines. Moti- 
vated by future instrumentation and recent proposals, recent 
studies have quantified the expected emission using cosmo- 



logical, hydrodynam ica l simulations jFurlanetto et aL 
iBertone et al]|2010al lbl: iTakei et al.ll201ll : iFrank et all 



2004 



2012 



We use a modifie d version of GADGET-3 (last described 
in ISpringej |2005| ). a smoothed particle hydrodynamics 
SPH) code that uses the entropy formulation of SPH 
Springel fc Hernauistll2002l ). which conserves both energy 
and entropy w here appropriate. This work is part of the 
OWLS project (jSchave et al.ll2010l ). which consists of a large 
number of cosmological simulations with varying (subgrid) 
physics. Note that the OWLS mo dels were also u s ed in t he 
studies of metal-line emission by IBertone et al] (|2010al lbh 
and IBertone fc Schavd (I2OI2D . Here we make use of the so- 
called 'reference' model, although we will compare the soft 
X-ray results with a simulation that includes AGN f eedback. 
The model is fully described in lSchave et al.l l|2010l ) and we 
will only summarize its main properties here. 

The simulations assume a ACDM cosmology with pa- 
rameters = 1 - 0.A = 0.238, flh = 0.0418, h = 0.73, 
as = 0.74, n = 0.951. These values are con sistent with the 
WMAP year 7 data (|Komatsu et al.ll201ll ). In the present 
work, we use the simulation output at redshifts 3 and 0.125. 

A cubic volume with periodic boundary conditions is 
defined, within which the mass is distributed over 512'' 



* The only significant discrepancy is in erg, which is 8 per cent, or 
2.3(T, lower than the value favoured by the WMAP 7-year data. 
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dark matter and as many gas particles. The box sizes 
(i.e. the length of a side of the simulation volume) of 
the simulations used in this work are 25, 50, and 100 co- 
moving /i~^Mpc. The (initial) particle masses for baryons 
and dark matter are 2.1 x 10*^ ( ^g h^iMpc )^ ^0 ^^'^ 1-0 ^ 
^^^^ 25 ft-°Mpc )^ respectively. The gravitational soften- 
ing length is 1.95 ( 25 fe-°Mpc ) comoving ft^^kpc, i.e. 1/25 of 
the mean dark matter particle separation, but is limited to 
a maximum value of 0.5 ( 25 h-°Mpc ) Proper h~^kpc, which 
is reached ai z — 2.91. 

Star formation is mo d elled according to the recipe 
of ISchave fc Dalla Vecchial l|2008l ). A polytropic equation 
of state Ptot oc pgas is imposed for densities exceeding 
riH = 0.1 cm~^, where Ptot is the total pressure and pg^s 
the density of the gas. Gas particles with proper densities 
riH ^ 0.1 cm~"^ and temperatures T ^ 10^ K are moved onto 
this equation of state and can be converted into star parti- 
cles. The star formation rate (SFR) per unit mass depends 
on the gas pre ssure and reprodu ces the observed Kennicutt- 
Schmidt law (jKennicuttl Il998l ) by construction. Particles 
with densities above nn = 0.1 cm~^ are not taken into ac- 
count for the results in this work. 

Feedback from star formation is impleme nted using the 
prescription o f iDalla Vecchia fc Schayj ((20081 ) . About 40 per 
cent of the energy released by type II supernovae is injected 
locally in kinetic form, while the rest of the energy is as- 
sumed to be lost radiatively. The initial wind velocity is 
600 km s~^ and the initial mass loading is two, meaning 
that, on average, a newly formed star particle kicks twice its 
own mass in neighbouring gas particles. 

The abundances of eleven elements released by mas- 
sive stars and intermediat e mass stars are followed as de- 
scribed in |w lersma et al.l (12009b'). Wc assume the stellar 
initial mass function (IMF) of Chabri er (2003) , ranging from 
0.1 to 100 Mq. As described in lWiersma et al] (|2009al '). ra- 
diative cooling and heating are computed element by ele- 
ment in the prese nce of the cosmic microw ave background 
radiation and the iHaardt fc Madaul (|200lh model for the 
UV/X-ray background from galaxies and quasars, assuming 
the gas to be optically thin and in (photo-)ionization equi- 
librium. 

Using the suite of OWLS models. Iwiersma et all (|201ll ') 
have shown that galactic winds driven by supernovae are es- 
sential for the enrichment of the IGM. Close to the centre of 
a halo, most gas has been enriched t o Z > 0.1 Zp), but at the 
virial radius the scatter is very large (|van de Voort fc Schavd 
|2012| ). 

Finally, we also show soft X-ray results for a simulation 
that includes AGN feedback (model AGN, Section [6]). Black 
holes grow via mergers and gas accretion and inject 1.5 per 
cent of the rest-mass energy of the accreted gas into the sur- 
rounding matter in the form of heat. The model is a mod ified 
version of the one introd uced by Springcl ct alj (|2005| ) and 
is described and tested in iBooth'fc Schavel (|2009l) . who also 
demonstrate that the simulation reproduces the observed 
mass density in black holes and the observed scaling rela- 
tions between black hole mass and both cent ral ste ll ar ve - 
locity dispersion and stellar mass. McCarth y et all (|201G| ) 
have shown that model AGN reproduces the observed stel- 
lar mass fractions, star formation rates, and stellar age dis- 



tributions in galaxy groups, as well as the thermodynamic 
properties of the intragroup medium. 

2.2 Identifying haloes 

The first step towards finding gravitationally bound struc- 
tures, is to identify dark matter haloes. These can be found 
using a Friends-of- Friends (FoF) algorithm. If the separation 
between two dark matter particles is less than 20 per cent 
of the average separation (the linking length 6 = 0.2), they 
are placed in the same group. Baryonic particles are linked 
to a FoF halo if their nearest dar k matter neighbour i s in 
that halo. We then use SUBFIND ('Springe l et al.l 120011 1 to 
find the most bound particle of a FoF halo, which serves as 
the halo centre and also corresponds to the location of the 
central galaxy. We compute the virial radius, 7?vir, within 
which the average density agrees with that predicted by 
the top-hat spherical collapse model in a ACDM cosmol- 
ogy (Brvan & Norman 1998). The halo mass, Mhaio, is the 
total mass contained inside Rvir- 

Table \T\ lists which simulations and which halo mass 
bins are used for various lines and redshifts. It also lists the 
median halo mass, median stellar mass, median star forma- 
tion rate, and the number of haloes in the different mass 
bins. 

2.3 Emission 

We calculate the emissivities of the gas following 
iBertone et al.l (|2010al ). We only summarize th e method here 
and refer the reader to iBertone et all (|2010al ') for details on 
the procedure. The names, wavelengths, and energies of the 
lines considered in this work are given in Table (2] 

We created emissivity tables as a function of tempera- 
ture, density, and redshift with the photo-ionization pack- 
age CLOUDY ( version c07.02.02), which was last described in 
iFerland et akl (fl998 V The gas is assumed to be optically thin 
and in ionization equilibrium i n the presence of t he co smic 
microwave background and the lHaardt fc Madaul (|200ll ) UV 
background. These assumptions were also made when cal- 
culating cooling rates during the simulation (|Wiersma et al.l 
i2009a ) and are thus fully consistent with the simulation. 
We further assumed solar abundances when creating the ta- 
bles, but the results are scaled to the abundances predicted 
by the simulation. The tables are created for temperatures 
10^ K< r < lO* '^ K with AlogioT = 0.05 bins and densities 
10"* cm"^ < wh < 10 cm"^ with Alogionn = 0.2 bins. 

We calculate the emission only for gas with nn < 
0.1 cm~'^, which is therefore by definition not star forming. 
Our simulations do not resolve the multiphase interstellar 
medium at nu 5^ 0.1 cm~'^. As we are interested in determin- 
ing the emission from diffuse gas in haloes, ignoring this gas 
is consistent with our aims. Note, however, that because of 
this omission, the mean surface brightness profiles presented 
in this work may underestimate the true mean profiles. 

The emissivity, e, for an emission line is derived from 
the tables as a function of logioT and logiown through log- 
linear interpolation. A gas particle's luminosity is then 

L = eiz,T,nn)^^ (1) 

in ergs~^, where mgas is its mass, p is its density, Xy is the 
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Table 1. List of the simulation box sizes and halo mass ranges used for different emission lines at different redshifts. Included are also 
the median halo mass, median stellar mass, median star formation rate, and the number of haloes in each halo mass bin. 



Redshift 


Band 


Simulation 


LogioMhaio range 


LogioMhaio 


LogioMstar 


SFR 


# haloes 








(Mo) 


(M©) 


(Mo) 


(Moyr-i) 




0.125 




REF 100 h^^Mpc 


12-13 


12.3 


10.8 


8.95 


2771 








13-14 


13.2 


11.8 


45.71 


301 








14-15 


14.2 


12.6 


144.25 


18 


0.125 


X-ray 


AGN 100 h-^Mpc 


12-13 


12.3 


10.3 


0.44 


2650 








13-14 


13.2 


11.1 


2.09 


242 








14-15 


14.2 


12.1 


48.67 


16 


0.125 


UV & Ha 


REF 50 h-iMpc 


11-12 


11.3 


9.3 


0.05 


2723 








12-13 


12.3 


10.9 


9.29 


374 








13-14 


13.2 


11.7 


38.58 


46 


3.0 


UV 


REF 25 /i-iMpc 


10-11 


10.3 


8.2 


0.11 


1857 








11-12 


11.3 


9.5 


4.24 


117 








12-13 


12.2 


11.0 


110.61 


5 



'SPH-smoothed' mass fraction of the element corresponding 
to the emission line, and X® is the solar mass fraction of 
the same el ement. We use SPH-sm oothed abundances as 
described bv lWiersma et al] (|2009bl ). This is consistent with 
the simulation, because these smoothed abundances were 
also used in the simulation to compute the cooling rates. 
The flux is 

(2) 



F = 



in erg s ^ cm ^ , or 



i-Kdf^ hpC 



+ 



(3) 



in photon s^^ cm~'^, with di^ the luminosity distance, hp the 
Planck constant, c the speed of light, and A the rest-frame 
wavelength of the emission line. 

To calculate the emission profiles, we project the flux 
using a flux-conserving SPH interpolation scheme onto a 
two-dimensional grid, centred on the halo centre. The sur- 
face brightness is calculated by dividing the flux by the solid 
angle subtended by a pixel either in arcsec^ or in sr, 

5b = F/n. (4) 

We will investigate the detectability of the emission 
lines listed in Tabled originating from the CGM. We chose 
these lines, because they are the brightest in the outer haloes 
of galaxies, although in the centres other lines may domi- 
nate. We will quantify the expected surface brightness as a 
function of radius for a range of halo masses. Our predictions 
for detectability should be considered with care. We did not 
take into account the velocities of the gas and thus the width 
of the lines and we ignored the effect of the Earth's atmo- 
sphere. Intervening absorption was also ignored. A robust 
study of detectability should produce datacubes, including 
noise, and analyze them using the same pipeline as the ob- 
servations. This is beyond the scope of the present work in 
which we show and discuss the theoretically expected sur- 
face brightness profiles. 

O VII is a line triplet at 0.574, 0.568, and 0.561 keV. The 
first line is a resonance line, the second one an intercombi- 
nation line and the last one a forbidden line. Photons at 



the resonance energy see larger optical depths and are thus 
more likely to be scattered in random directions, giving rise 
to lower observed fiuxes or even absorption. The intercom- 
bination line is significantly weaker (0.5 — 0.6 dex) than the 
other two. The resonance line is only slightly stronger than 
the forbidden line (by about 0.1 dex) without taking into 
account attenuation. Although the actual 0.574 keV surface 
brightness can in principle be much lower due to resonant 
scattering, we estimated the maximum optical depth for the 
O VII resonance line in a halo to be very small. Out of the 
haloes we considered, only 20 have at least one pixel with 
optical depth r > 0.1 (using Equation 44 in iKaastra et al.l 
[2008 with a velocity dispersion of 100 km s"'^). The line is 
optically thin in the vast majority of cases, so we do not 
expect resonant scattering to infiuence our results. Another 
effect that could play a role is charge exchange, caused by 
collisions between hydrogen atoms and oxygen ions, in which 
case one would find enhanced emissi on for the forbidden 
O VII line (e.g. iLiu et all [20T1I . I2OI2I ). We only show the 
surface brightness of the forbidden line at 0.561 keV, but 
note that if the emission at all three wavelengths is added 
together or if charge exchange is important, then the signal 
could be more than twice as strong. 

The UV doublets, C iv, O vi, and Si iv have flux ratios 
of 2:1. We consider only the strongest line in this paper, but 
the results for the weaker lines thus follow directly. 



3 EXPECTED DETECTION LIMITS 
3.1 Soft X-ray 

Current X-ray telescopes, i.e. Chandr(^ XMM-Newtot^ 
and Suzaki^ have already detected continuum and metal- 
line emission from the diffuse, intracluster gas. To map the 
halo gas out to larger radii and in lower mass haloes with 



http: / /chandra. harvard.edu/ 



^ http://xmm.esac.esa.int/ 
■* http:/ /www .jaxa.jp/projects/sat/astro_e2/index_e. html 
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Table 2. List of ion, rest-frame wavelength, and rest-frame energy 
for the soft X-ray, optical, and UV emission lines used in this 
work. 



ion 


A (A) 


E (eV) 


C VI 


33.74 


367 


N vii 


24.78 


500 


O VII 


22.10 


561 


O VIII 


18.97 


654 


Ne X 


12.14 


1021 


H I (Ho) 


6563 


1.89 


C III 


977 


12.69 


C IV 


1548 


8.01 


O VI 


1032 


12.01 


Si III 


1207 


10.27 


Si IV 


1394 


8.89 



future instruments, a combination of high angular and spec- 
tral resolution and a large field of view are preferred (e.g. 
iBertone et al.l[2010al ). Higher spectral resolution results in 
higher signal-to-noise for unresolved lines. As the emission 
is dominated by bright, compact regions, high angular res- 
olution is necessary to prevent smearing out the emission, 
in which case it would be harder to detect as well as im- 
possible to correctly identify its origin. ASTRO-t^ will be 
launched in a few years. While its spectral resolution is rel- 
atively good (5 eV), its field of view is rather small (3x3 
arcm in^) and its spati al resolution is coarse (1.3 arcmin) 
(e.g. ISoong et al.l l201ll . see also Table (H)). Therefore, it is 
not an ideal instrument for studying halo X-ray profiles, 
although it will do relatively well for large targets, such as 
massive clusters. A number of X-ray missions have been pro- 
posed inJhe_2as^J[eaxSj^uch as IXO/ATHEN^ and ORI- 
GIN l|den Herder et aLll201ll '). Their proposed specifications 
are listed in Table |31 as are those for existing and future UV 
and optical instruments, which are discussed in Sections l5.2l 
and 15.31 respectively. 

For CRIS on ORIGIN, the instrumental background is 
2 X lO"'^ photon s~^cm~'^keV~^, which corresponds to 1.3 x 
10^ photon s~^sr~^keV~^ (with a pixel size of 300 pim and 
24 arcsec) . In the narrow-field mode, the energy resolution is 
2.5 eV, so the instrumental background is 3 photon s~^sr~^. 
For an effective area of 10"^ cm^, the surface brightness of 
the instrumental background thus equals 3 x 10~^ photon 
s~^cm~^sr~^. (It would be twice as high in the wide-field 
mode.) 

For XMS on IXO/ATHENA the instrumental back- 
ground is also 2 X 10~^ photon s~^cm~^keV~^, which cor- 
responds to 8.5 X 10* photon s~^sr~^keV~^ (with a pixel 
size of 300 /.im and 3 arcsec). For an energy resolution of 
2.5 eV, which is appropriate for the narrow-field mode, this 
corresponds to 2 x 10^ photon s~^sr~^. For an effective area 
of 3 X 10* cm^, the surface brightness of the instrumental 
background thus equals 7 x 10~^ photon s~^cm~^sr~^. (It 
would be four times as high in the wide- field mode.) 

For SXS on ASTRO-H the instrumental background 

^ http://astro-h.isas.jaxa.jp/ 

^ |http:/ /ixo.gsfc.nasa.govTI http:/ /sci.esa.int/ixo 



is an order of magnitude lower, 2 x 10^'^ photon 
s^^cm~^keV~^, which corresponds to 7.0 x 10^ photon 
s~^sr~^keV~^ (with a pixel size of 830 /xm and 29 arc- 
sec). For an energy resolution of 5 eV, this corresponds 
to 3 photon s~^sr~^. For an effective area of 50 cm^, the 
surface brightness of the instrumental background equals 
7x 10~^ photon s~'^cm~^sr~^, an order of magnitude higher 
than for XMS on IXO/ATHENA. 

For XIS on Suzaku the instrumen tal background is 
higher, ~ 2 x 10 photon s~^cm~^keV"^ (|Tawa et al.ll2008l ). 
but because it contains many pixels, this corresponds to 
only 4.2 x 10^ photon s~^sr~^keV~^ (with a pixel size of 
24 fim and 1 arcsec). Given an energy resolution of 70 eV, 
this amounts to 30 photon s~^sr~^. For an effective area 
of 100 cm^, the surface brightness of the instrumental back- 
ground equals 3 x 10~^ photon s^^cm^'^sr"*, a factor of four 
higher than for SXS on ASTRO-H. 

iHickox fc MarkevitchI (|2006l ) determined that the un- 
resolved cosmic X-ray background is about 14 ± 1 photon 
s~^cm~^sr~^keV~* in the energy range we are interested 
in (0.5 — 1 keV). This unresolved signal is thought to be 
dominated by diffuse Galactic emission and local thermal- 
like emission. Part of this unresolv ed background could 
be re solved in deeper observations l|Hickox fc MarkevitchI 
I2OO7I ). We will conservatively use the limit of 14 photon 
s~*cm^^sr~^keV^^ as it is an upper limit on the unresolved 
background for future observations. For the narrow-field 
mode, the energy resolution of CRIS on ORIGIN and XMS 
on IXO/ATHENA (i.e. 2.5 eV) this corresponds to about 
3x 10~^ photon s~^cm~^sr~'^. The instrumental background 
is therefore lower than the unresolved X-ray background. 
Observations should be able to probe down to 10" * pho- 
ton s~^cm~'^sr~^ without being strongly affected by the 
background. For this surface brightness a 1 Msec expo- 
sure would result in about 2 detected photons per resolu- 
tion element with ORIGIN at 30 arcsec resolution and with 
IXO/ATHENA at 5 arcsec resolution. 

The energy resolutions of SXS on ASTRO-H and of the 
wide-field mode for CRIS on ORIGIN are two times worse 
than the narrow-field mode for CRIS and XMS. The unre- 
solved cosmic X-ray background is therefore 7 x lO"'^ pho- 
ton s~*cm~^sr~^, comparable to the instrumental back- 
ground of SXS on ASTRO-H. The wide-field mode for XMS 
on IXO/ATHENA leads to unresolved cosmic X-ray back- 
grounds of 14 X 10"^ photon s"^cm~^sr~^. Observations 
can therefore only probe down to several times 10~^ pho- 
ton s~^cm~'^sr"^ without being strongly affected by the 
background. Assuming a detection limit of 3 x 10~* pho- 
ton s~'^cm~^sr~'^, a 1 Msec exposure would result in about 
2 detected photons per resolution element on ASTRO-H at 
78 arcsec resolution. 

The energy resolution of XIS on Suzaku is only about 70 
eV, as this is an imager and not a spectrograph. Therefore, 
the unresolved X-ray background corresponds to 1 photon 
s~^cm~^sr~^, a factor of three larger than its instrumental 
background. Observations can only probe down to about 
5 photon s~^cm~'^sr~^ without being strongly affected by 
the background. Additionally, the low energy resolution, 70 
eV, may lead to strong projection effects and would make 
identifying spectral lines very difficult. 
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Table 3. Summary of the most relevant technical specifications for instruments discussed in the main text. 



Telescope (instrument) 


Field of view 


Angular resolution 


Spectral resolution 


Effective area 


Instrumental background 




(arcmin^) 


(arcscc) 




(cm^) 


(s-icm-^keV-i) 


ORIGIN (CRIS) 


10x10 (30 X 30) 


30 


2.5 (5) eV 


10^ 


0.02 


IXO/ATHENA (XMS) 


2x2 (5 X 5) 


5 


2.5 (10) eV 


3 X 10* 


0.02 


ASTRO-H (SXS) 


3x3 


78 


5 eV 


50 


0.002 


Suzaku (XIS) 


18 X 18 


120 


70 eV 


100 


2 X 10 


FIREBALL 


2.67 X 2.67 


10 


0.4 A 


8 X lO"* 




Hale (CWI) 


1 X 0.67 


2.5 X 1 


0.12 A 


2 X 10^ 




Keck (KCWI) 


0.33 X (0.14 - 0.57) 


0.35 - 1.4 


0.03-0.6 A 


8 X 10^ 




VLT (MUSE) 


1 X 1 


0.2 


0.4 A 


5 X 10^ 





For reference, 

in-19 -1 -2 -2 

10 ergs cm arcsec 

.2 / A \ / 1 + z 



2.7 X 10' 



103 A 



Y2^^ photons ^cm ^sr ^. (5) 



3.2 Low-redshift UV 



To detect UV metal-line emission as a tracer of halo gas, in- 
struments should ideally have a large field of view and a high 
spatial a nd angular resolution, as is also the case for X-ray 
emission l|Bertone et al.|[2010bl ). The emission will be domi- 
nated by relatively high-density material. With high spatial 
resolution, it will be possible to measure its dumpiness. 

A surface brightness limit of order 10~^* erg 
s~^cm~^arcsec~^ , as is appropriate for th e FIREBALj^ bal- 
loon experiment (|Tuttle et al.ll2008l . [201ol ). will only be suffi- 
cient to detect metal-line emission in the centres of massive 
haloes. The wavelength range of FIREBALL (2000-2200 A) 
is such that it will only probe the UV metal lines at a some- 
what higher redshift than shown in this paper (z ~ 0.35 for 
C iv) where the surface brightness is a bit lower. A detection 
limit of order 10"'^^ erg s~^cm~^arcsec~^ is en visioned for 
the next generation UV mission ATLAST (|Postman.200d ). 
which would be a huge improvement. 

For reference. 



10-19 -1 -2 -2 
10 ergs cm arcsec ~ 

8.5xl0^f^Vi±i) ph, 
Vi03A/V 1.25 / ^ 



lotons ^cm ^sr ^. (6) 



3.3 High-redshift UV 

It is possible to detect rest-frame UV emitted at high red- 
shift from the ground, in the optical. In the near future, 
several integral field unit spectrographs will become op- 
erationa l, such as the Mu lti-Unit Spectroscopic Explorei[f| 
(MUSE jSacon et al.ll2O10l ) on the Very Large Telescope and 
the Keck Cosmic Web ImageiQ (KCWI; iMartin et aLll2010l) 
on Keck. The Cosmic Web Imagei[3 (CWI; Rahman et al.l 



http://www.srl.caltech.edu/sal/fireball.html 
* http://muse.univ-lyonl.fr/ 

^ http://www.srl.caltech.edu/sal/keck-cosmic-web-imager.html 
http: / / |htt p : / / www . srl . caltech . edu/sal/ cosmic- web-imager . html | 



I2OO6I : iMatuszewski et al.1 [2OI0I ) is already installed on the 
200-inch Hale telescope. 

With narrow-band observations it has already been 
demonstrated that diffuse H I Lya emission can be detected 
down to ~ 10~^* erg s~ ^ cm~^ arcsec"^ in indi vidual objects 



(e.g. ISteidel et al]|200Gl : iMatsuda et"al]l2004l) and down t o 
~ ]^g-i9 gj.g s~^cm~'^arcsec~^ in stacks ( Steidel et al ] |201ll ). 
CWI is expected to reach similar depths. New instrumenta- 
tion, such as MUSE and KCWI, will allow us to go one order 
of magnitude deeper. This incre ase in sensitivity is essen tial 
for the detection of metal lines (|Bertone fc Schav3l2012l ). 



4 SURFACE BRIGHTNESS PROFILES 

For all the emission lines listed in Tabled we calculate sur- 
face brightness maps for each halo and stack the maps for 
haloes within a certain halo mass range. We then average the 
mean surface brightness maps azimuthally to make emission 
profiles. In this Section we discuss the surface brightness 
profiles and the detectability of the different emission lines. 
Median profiles are calculated by taking the median of the 
pixels in the stacked images and then averaging those az- 
imuthally. 



4.1 Soft X-ray 

In Figure [T] we show the mean (median) X-ray line surface 
brightness profiles for haloes at z = 0.125 as solid (dotted) 
curves. We use three halo mass bins, roughly corresponding 
to (the haloes of) galaxies, groups, and clusters, with masses 
of Mhaio = 10'^-'^ Mq, 10"-'* M©, and lO'"*-'"^ M©. We 
show the profiles out to twice the median virial radius and 
with a thickness of four times the maximum virial radius in 
each mass bin. A spectral resolution of 2.5 eV corresponds 
to a thickness of about 15 Mpc. Using 15 Mpc thick slices 
would only mildly fiatten the profiles of the lowest mass 
haloes around i?vir, which is below the detection limit of 
proposed instruments. 

The strongest observable soft X-ray line is O vill (0.654 
keV). It is followed by C vi (0.367 keV) and O vil (0.561 
keV) for galaxies. Nex (1.021 keV) has the lowest surface 
brightness for A/haio ~ 10^^~^^ M©. In groups, Ne x is the 
second strongest line in the centre of the halo, but it is 
weaker than C vi and O vil in the halo outskirts. As men- 
tioned in Section [2. 3 1 O vil is a triplet and the emission will 
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Figure 1. Mean (solid curves) and median (dotted curves) surface briglitness as a function of radius at 2: = 0.125 for tlie soft X-ray 
lines indicated in the legend. We use three different halo mass bins with masses of M^^^i^ = 10^^~^^ Mq (left), 10^'^"^'' Mq (middle), 
and 10^^~^^ Mq (right). The thickness of the slices is 1.7, 3.7, and 5.6 comoving /i~^Mpc, respectively, corresponding to four times the 
largest virial radius in each mass bin. The radius on the bottom x-axis is normalized by dividing by the median virial radius in each 
mass bin. The surface brightness in photon s^^cm^^sr"'^ indicated by the left y-axis is exact. For the right y-axis it has been converted 
to erg s^^cm^'^arcsec"^ using (A) = 22.2 A (0.558 keV). The pixel size was 5 arcsec (8.9 comoving h~^kpc; 7.9 proper /i~^kpc) before 
binning. Ovill is the brightest line at all halo masses, followed by Cvi. The relative strengths of the lines vary with halo mass. The 
profiles flatten at _R < 10 /i~^kpc, because this region is dominated by the ISM, which we excluded from our analysis. 



be about a factor of two stronger when the contributions of 
the three lines are added. For both galaxies and groups, the 
Ne X profile is steeper than that of all the other lines shown. 
In clusters C vi and Ne x are of comparable brightness at 
aU radii and about 0.6 dex below O vill. N vil (0.500 keV) 
is similar to O vil at all radii for cluster haloes and about an 
order of magnitude fainter than O vill. We also computed 
the emissivities of C v (0.308 keV), N vi (0.420 keV), Ne ix 
(0.922 keV), Mg xii (1.472 keV), Si xiii (1.865 keV), Si xv 
(2.460 keV), and Fe xvii (0.727 keV), which were in genera l 
weaker than the lines shown (see also lBertone et al]|2010al V 
The surface brightnesses of Ne ix and Fe xvil are as strong 
as those of N x in galaxies and Fe xvil is as high as O vil 
for R < 0.2i?vir in groups. Both lines are fainter in cluster 
centres, but possibly still detectable. 

The median surface brightness profiles are lower than 
the mean profiles at all radii. A large difference between 
mean and median profiles indicates that the emission is 
dominated by relatively rare, bright pixels and thus that 
the emission is clumpy. Note that the median profiles have 
still been averaged azimuthally. For Mhaio ~ 10^^"^"^ Mq, 
Ne X and O vill exhibit large differences between mean and 
median profile (the median Ne x profile is lower than the 
range shown), whereas for Mhaio = 10^*~^^ Mq this is the 
case for O vil. The difference increases at very small radii 
{R < 0.1i?vir) for Mhaio = 10^^"^^ Mq and towards the edge 
of the halo {R > 0.57?vir) for all haloes. At intermediate radii 
the difference between the mean and median profiles is small 
for the majority of the lines, implying that the X-ray emis- 
sion is relatively smooth. 

At a limiting surface brightness of 10~^ photon 
s^^cm^^sr^^, as is appropriate for CRIS on ORIGIN and 
XMS on IXO/ATHENA, we would be able to detect O vill 
emission out to 80 per cent of the virial radius of a clus- 
ter with Mhaio = 10""^^ Mq, C VI out to 40 per cent, 
and O VII and N vil out to 20 per cent. For groups, 
Mhaio ~ 10^^~^* Mq, O vii can be observed out to the same 



physical scale, ~200 comoving h^^kpc, which corresponds 
to a larger fraction of the virial radius, 0.47?vir. In this case, 
O vill is observable out to 0.7i?vir and C VI out to 0.5i?vir. 
For Mhaio = 10^^"" Mq, O viii, O vii, and C vi can be 
observed out to the same radius of 0.3i?vir. 

SXS on ASTRO-H has poor angular resolution (78 arc- 
sec which equals 0.14 comoving /i~^Mpc), but it may be 
able to detect Ovill emission in the centres of groups and 
clusters. It will not be able to see detailed structure, because 
the emission would be contained in a single pixel. 



4.2 Low-redshift UV 

Figure [5] shows the surface brightness for the UV metal lines 
listed in Table [2] at 2 = 0.25 for haloes with masses of 
Mhaio = 10"-i^ Mq, 10^^-^^ Mq, and 10'^-''' Mq. Line 
styles are identical to those in Figure [T] (i.e. solid indicates 
means, dotted medians), but colours show different emission 
lines, as indicated in the legend. We also include hydrogen 
Ha (black, dashed curves) from the same haloes and discuss 
it in Section WM 

C in (977 A) is the brightest line, followed by Si in 
(1207 A). The emission profile for all lines peaks in the 
halo core, where it may, however, be outshone by the cen- 
tral galaxy (recall that we omitted emission from gas with 
uh > 0.1 cm~^) O VI (1032 A) gives the most extended pro- 
file, in the sense that the ratio between emission in the halo 
core and outskirts is smallest. This is caused by the fact that 
O VI traces hotter and more diffuse gas than the other lines. 
C IV (1548 A) and Si iv (1294 A) are weaker than C III and 
Si in, but because they have frequencies redward of Lya, 
they are less easily absorbed by intervening gas clouds. We 
also computed the emissivities of He 11 (1640 A) and N v 
(123 9 A), which w ere we aker and therefore not shown (see 
also iBertone et aU bOlObl ) . 

At all radii, the median surface brightness profiles are 
lower than the mean profiles. For Mhaio = 10^^~^^ Mq, 
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Figure 2. Mean (solid curves) and median (dotted curves) surface briglitness profiles for tiie emission lines and halo masses indicated 
in the legend at z = 0.25. Note that the halo masses, M^aio = lO^^-^^ Mq (left), IQl^-lS M0 (middle), and lO^^-^* Mq (right), are 
lower than in Figure[T] The dashed curve shows the median surface brightness profile for Ha. The thickness of the slices is 0.8, 1.7, and 
3.3 comoving /t~^Mpc, respectively, corresponding to four times the largest virial radius in each mass bin. The surface brightness in 
erg s~^cm~'^arcsec~^ on the left ly-axis is exact. For the right y-nxis, it has been converted to photon s~^cm~^sr~^ using (A) = 1232 A. 
The pixel size was 2 arcsec (3.5 comoving h~^kpc] 2.8 proper /i~^kpc) before binning. Cm is the brightest UV metal line, followed by 
Sim. 



all lines show a large difference between mean and me- 
dian profile, especially in the inner parts, by up to 1.5 dex. 
For Mhaio = 10^^"^^ M0 the difference is very small at 
0.057?vir < R < 0.4i?vir- C III, Si III, and Si iv have similar 
mean and median surface brightness profiles at ii < 0.1i?vir 
for Mhaio = iQ^^-^'i Mq, whcrcas for O vi this is the case 
at _R > 0.3i?vir. The mean and median C iv profiles are dif- 
fer by about 0.5 dex. For Mhaio > 10^^ M©, O vi exhibits 
the smallest difference between mean and median profiles, 
consistent with the expectation that it traces hotter, more 
diffuse gas. 

A surface brightness limit of 10~^* erg s~^cm~'^ arcsec"^ 
should enable us to observe metal-line emission from the cen- 
tres of fairly massive galaxies (Mhaio > 10^^ M0. However, 
it is possible that the central galaxy will outshine the halo 
in the halo cores. 

By stacking haloes, FIREBALL could in principle 
probe down to a surface brightness limit of 10~^^ erg 
s~^cm~^arcsec~^. Its wavelength rage is such that it can 
only detect Cm and Sim at 2: ~ 1. As the surface bright- 
ness is ~ 0.3 dex lower at 2: = 1 as compared to 2: = 0.25 
(not shown), it may see emission out to ~ O.liivir for 
Mhaio > 10^^ Mq . A complicating factor is that C ill, O vi 
and Sim have frequencies blueward of Lya and their light 
could therefore be absorbed by intervening hydrogen, an ef- 
fect that is largest for Cm and smallest for Sim. Civ and 
Si IV have frequencies redward of Lya and are thus not af- 
fected by this. 

With a detection limit of order 10~^^ erg 
s~^cm~'^arcsec~^ we would be able to detect Cm out 
to the virial radius and several other metal lines out to 50 
per cent of i?vir- Such a low surface brightness limit may be 
attainable for t he next generation UV mission, ATLAST 
jPostmanl |2009| V As Lya emission is even brighter, such 
a mission can even probe the gas in emission outside 
galactic haloes, in th e filaments of the cosmic web (e.g. 
iFurlanetto et al.ll2003l '). 



4.3 High-redshift UV 

Figure [3] shows the UV line surface brightness profiles at 
2 = 3 for haloes with masses of Mhaio ~ 10^""^^ Mq, 
10""^^ Mq, and 10^^"^^ Mq. The surface brightness is 
about an order of magnitude lower at 2 = 3 than at 2 = 0.25 
at a fixed fraction of -Rvir- Additionally, the profile shows 
a strong fiattening inside O.l-Rvir- This is because we are 
probing down to smaller physical scales, where the emission 
is dominated by star-forming gas, which we have excluded 
from our calculations. 

Comp aring our metal-line results to the Lya profile ob- 
served by ISteidel erahl (|201ll ) for Mhaio ~ 10^^ Mq and 
(2) = 2.65, we predict that the C lll surface brightness is 
about 1-2 dex lower than the Lya brightness (ignoring ab- 
sorption by intervening gas). The other metal lines that are 
shown are about 0.5-1 dex fainter than C ill. Because C ill 
has a much shorter wavelength than Lya, it is strongly ab- 
sorbed by the intervening medium. It is therefore possible 
that the other metal lines will be observed to be stronger 
than C III. We calculated, but do not show, the Lya surface 
brightness profile from our simulations in the optically thin 
limit and without a contribution from star formation. The 
observed profile is ~ 0.7 dex higher than the predicted pro- 
file, but this is not surprising considering that most of the 
Lya emission is thought to be originating from H 11 regions 
in the ISM (e.g. IFurlanetto et al.ll2005l : ISteidel et al]|201ll ') 
and that we excluded this gas from our analysis. 

For Mhaio = 10^°"" Mq, all lines show a 0.3-1 dex 
difference between mean and median profiles. The differ- 
ence is largest in the inner parts. For Mhaio > lO'^^ Mq 
at intermediate radii (0.03-Rvir < R < O.Gflvir) this differ- 
ence is maximally 0.3 dex. The difference increases towards 
the edge of the halo. The UV emission is least clumpy for 

Afhalo > 10^' Mq. 

A detection limit of 10~^^ erg s~^cm~'^ arcsec"^, which 
should be possible for CWI, is not low enough for the de- 
tection of most metal lines. Cm emission could be seen 
out to almost 0.1 /ivir for massive haloes, if it is not at- 
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Figure 3. Mean (solid curves) and median (dotted curves) surface brigtitness profiles for the rest-frame UV emission lines and halo 
masses indicated in the legend at z = 3. Note that the halo masses, M^^^i^ = 10^'^'^^ Mq (left), lO^^'^^ Mq (middle), and 10^^-^^ Mq 
(right), are lower than in Figures [TlandlSI The thickness of the slices is 0.45, 0.95, and 1.25 comoving h^'^Mpc, respectively, corresponding 
to four times the largest virial radius in each mass bin. The pixel size was 0.2 arcsec (4.6 comoving h^'^kpc; 1.2 proper h~^\cpc) before 
binning. Cm is the brightest line. Civ, Sim, and Si IV are of similar strength. The surface brightnesses are lower than at z = 0.25 and 
the flattening at small radii is more pronounced. 



tenuated by the IGM. With a detection limit of 10~ erg 
s~^cm~'^arcsec~^, which should be possible for MUSE and 
KCWI when stacking many galaxies, C ill can be seen out 
to 0.3i?vir for Mhaio = 10""^^ Mq and out to 0.2i?vir for 
-|^q12-13 Mq at z = 3. Note, however, that MUSE can 
only see C ill emission for z > 3.8. Other metal lines are 
weaker, but can still be probed to 10 per cent of J?vir for 
Mhaio > lO'^^ Mq . The emission may, however, in reality be 
brighter at these small radii, because we excluded emission 
from gas with nn > 0.1 cm"''. We predict that the sur- 
face brightness for lower-mass haloes is below this detection 
limit. 



4.4 Low-redshift Ha 

Hq (6563 A) emission at low redshift is observable with the 
same instruments as high redshift rest-frame UV emission. 
The predicted emission from gas around galaxies is shown 
in Figure [2] as a dashed, black curve. Interestingly, the Ha 
emission is comparable to C lll emission at 7? > 50 kpc for all 
haloes with Mhaio = lO""" Mq. However, for i? < 50 kpc, 
Hq is much brighter than C ill for the lowest halo mass bin 
(Afhaio = 10^^~^^ Mq) and up to an order of magnitude 
fainter for the highest halo mass bin (Mhaio = 10"-" Mq). 
We stress again that we ignored emission from gas with tih > 
0.1 cm"^, implying that the true emission may be brighter. 

For Mhaio ~ lO^^-^'^ Mq the median surface brightness 
profile is significantly lower than the mean profile (2 dex at 
0.2_Rvir), but this difference decreases towards smaller radii. 
This indicates that the emission is more clumpy in the outer 
than in the inner halo. For higher-mass haloes, the difference 
is smaller, but also increases at large radii. 

Assuming a limiting surface brightness of 10"^*^ erg 
s~^cm~^arcsec~^. Ha emission from halo gas can be de- 
tected out to 0.5iivir for haloes with Mhaio ~ 10^^"^^ Mq, 
to 0.25J?vir for haloes with Mhaio = 10^^"" Mq, and out 
to 0.47?vir for group-sized haloes. This may be feasible when 
stacking deep observations centred on galaxies with MUSE 
or KCWI. 



5 FLUX- WEIGHTED PHYSICAL 
PROPERTIES 

To investigate possible observational biases when deriving 
physical properties of the halo gas. Figures |4}{6] show the 
logarithms of the flux-weighted median temperature, over- 
density, and metallicity as coloured curves for the same red- 
shifts and halo mass bins as in Figures[T]-[3] The legends indi- 
cate the difi'erent emission lines. Mass- and volume- weighted 
properties are shown as black and grey curves, respectively. 
The dotted curves show the 16th and 84th percentiles and 
the solid curves show medians. The top panels show the 
results &t R < 0.5i?vir, while the bottom panels show 
0.5i?vir — -Rvir- These radii are the true three-dimensional 
radii, but the results change by less than 0.1 dex if we use 
projected radii (not shown). 

5.1 Soft X-ray 

Figure |4] shows fiux-weighted properties at z = 0.125 for soft 
X-ray lines. The temperature at which the emissivity peaks 
is highest for Ne x (« lO*^'* K) and O vill (« 10*^'^ K), and it 
de creases slowly toward s higher temperatures (e.g. Figure 1 
of iBertone et al. | |2010al ). This is the reason why, for galax- 
ies and groups, these lines have the highest flux-weighted 
temperatures. For clusters, the Ne x-weighted temperature 
is lower than that of O vill and close to the temperature 
where its emissivity peaks. The fiux-weighted median tem- 
peratures increase less steeply with halo mass than the mass- 
weighted median temperature, although they still increase 
by ^ 0.2 — 0.9 dex for an increase of 2 dex in halo mass. 
For the lowest (highest) mass bin the fiux-weighted temper- 
atures are higher (lower) than the mass-weighted temper- 
ature. The Ovil fiux-weighted temperature stays roughly 
constant with halo mass. This is due to the fact that the 
emissivity curve of O vil is more strongly peaked in tem- 
perature th an the e missivity cur ves of the other lines (see 
Figure 1 of iBertone ct al. 2010i). It peaks at lO*^'^ K and 
drops quickly for both higher and lower temperatures. 

In the halo outskirts, the temperatures probed are 
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Figure 4. Mass- weighted (blaclt curves), volume weiglited (grey curves) and fiux-weigiited (coloured curves for different soft X-ray 
emission lines as indicated in the legend) median temperature, overdensity, and metallicity at 2 = 0.125 as a function of halo mass. 
The top panels show the results for gas at i? < 0.5-??vir> while the bottom panels show 0.5-Rvir — Rviv The dotted curves show the 
mass- and volume-weighted 16th and 84th percentiles. The flux-weighted properties are biased towards high density and metallicity. The 
flux-weighted temperature is biased towards the temperature at which the line reaches its peak emissivity. 



somewhat lower, by 0.2 — 0.3 dex, than in the inner parts. 
This is true for flux- weighted as well as mass- weighted me- 
dian temperatures and densities, so they follow the same 
trend. 

The densities are 0.5 — 0.7 dex lower in the outer halo 
than in the inner halo. The flux- weighted overdensities of 
the gas responsible for the different lines are very similar 
and close to the 84th percentile of the mass-weighted mean 
overdensity which is ~ 0.2 — 0.7 dex above the median mass- 
weighted density. The reason for this is that the emissivity 
scales as p^, which biases the emission towards high densi- 
ties. The flux-weighted densities increase by ~ 0.8 dex over 
two orders of magnitude in halo mass, but the volume- and 
mass- weighted median densities increase by only 0.2-0.5 dex 
over the same range of halo masses. 

The median mass-weighted metallicity is 0.04-0.2^0, 
but the flux- weighted metallicities are about 0.5-1.0 dex 
higher. For galaxies and groups, the median flux-weighted 
metallicities are similar to the 84th percentile of the mass- 
weighted metallicity. The flux-weighted metallicities for clus- 
ters are signiflcantly above the 84th percentile, especially for 
O VII, which also has the highest overdensity and lowest tem- 
perature. Thus, not surprisingly. X-ray metal-line emission 
is biased towards high-metallicity gas. As metal lines domi- 
nate the emission in the soft X-ray band (0.5-1. kcV), the 
same w ill be true for broad band emission (e.g. ICrain et al.l 
l2010al ). 



For clusters, we also flnd significant bias for line-of-sight 
velocities and velocity dispersions as measured with the soft 
X-ray lines described above (not shown). The velocity bias 
is highest for O vil, which is also most discrepant in other 
properties. For harder X-ray lines, such as Fe xxv, whose 
emissivity curves peak at higher temperatures, the bias for 
all quantities decreases and almost disappears for the flux- 
weighted velocities and velocity dispersions. We therefore 
conclude that it is essential to use metal lines with peak 
emissivity temperatures close to the virial temperature of 
the halo when deriving gas properties. We caution, however, 
that he bias does not vanish completely, so it is important 
to compare properties derived from real and simulated ob- 
servations. 



5.2 Low-redshift UV 

Figure [5] shows flux-weighted properties at 2 = 0.25 
for UV lines. The flux-weighted temperature depends on 
the speciflc metal line. They are insensitive to the mass- 
weighted temperature of the halo, because t he emissivities 
are str ongly peaked at specific temperatures (|Bertone et al.l 
l2010bl ). Hence, detection of rest-frame UV metal-line emis- 
sion puts strong constraints on the temperature of the emit- 
ting gas. The difference between Sim and Ovi is just over 
an order of magnitude. All UV metal-lines probe gas with 
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temperatures lO*'^ * K with the exception of O vi, which is 
dominated by 10*'^ K gas. 

The different metal hnes probe very different densities, 
all higher than the mass-weighted overdensity and increas- 
ing as the corresponding temperature decreases. O VI results 
in the lowest flux-weighted overdensity. The flux-weighted 
density increases via C IV to C ill to Si iv and is highest 
for Sim. For all halo masses, the silicon lines have higher 
flux- weighted densities than the 84th percentile of the mass- 
weighted distribution, fn all haloes the flux-weighted median 
overdensities are higher than the volume- weighted 84th per- 
centile. In group-sized haloes all lines are dominated by gas 
even denser than the mass-weighted 84th percentile. The 
difference between the O vi and Si ill flux- weighted overden- 
sities is 1 dex. The flux-weighted overdensities increase 
more steeply with halo mass than the volume- and mass- 
weighted overdensity, by about f.5 dex for an increase of 
2 dex in halo mass. The overdensity is consistently lower by 
0.3 — 0.8 dex in the halo outskirts. 

The flux-weighted median metallicity is 0.5-f.5 dex 
higher than the mass-weighted median metallicity and it is 
also higher than the mass-weighted 84th percentile. The dis- 
crepancy is larger in the halo outskirts. For high-mass haloes 
with A/halo > fO^^ M0, O VI probes gas with ~ 0.3 dex 
lower metallicities than the other lines. In the inner parts of 
Afhaio > 10^^ Mq haloes, the emission of all lines is domi- 
nated by gas with supersolar metallicities, even though on 
average the gas is only enriched to « O.f Zq. All metal- 
licities increase with halo mass. The metallicity difference 



between inner and outer halo is largest for the lowest-mass 
haloes. 



5.3 High-redshift UV 

In Figure [5] the flux- weighted properties are shown at z = 3 
for UV lines. The general picture is unchanged at high red- 
shift as compared to low redshift (see Figure [5]). The flux- 
weighted temperatures are insensitive to halo mass. They 
trace temperatures near or above the median mass- weighted 
temperature for low- mass haloes (Mhaio < 10^^ M© ), but 
the opposite is true for higher-mass haloes. The scatter in 
temperature increases with halo mass and becomes very 
large, with the mass-weighted 84th percentile 2.5 dex above 
the f6th percentile. This reflects the bimodal nature of 
halo gas in massive high-redshift galaxies, where the gas 
either has a temperat ure close to the virial tempera ture or 
has T < 10^ K (e.g. Ivan de Voort fc Schavd [20]! ) . O vi 
traces gas about an order of magnitude warmer than the 
other metal lines shown, because its emissivity peaks around 
10^'^ K compar ed with 10"^-^"^ ^.fa ^j ^*^ other metal hnes 
(see Figure f of lBertone fc SchavdlioH ). 

The density differences between different lines is some- 
what reduced at high redshift as compared to low redshift. 
The flux-weighted overdensities are lower, but the corre- 
sponding physical densities are higher at high redshift. 

At z = 3 the scatter in the density is larger than at 
z — 0.25: the difference between the I6th and 84th percentile 
is more than an order of magnitude as opposed to 0.4 dex 
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at z = 0.25. As a result, none of the lines are dominated by 
overdensities above the 84th percentile at high redshift. The 
flux-weighted densities increase with halo mass, even if the 
volume- and mass- weighted densities do not. The median 
densities are lower in the outskirts than in the centre, by up 
to ~ 0.6 dex. 

As at low redshifts, the flux-weighted properties are 
biased towards high metallicities. All lines trace gas with 
metallicities of about 0.1-0.5 Zq. The mass- weighted metal- 
licity is up to 1.2 dex lower in the outer halo than in the inner 
halo, whereas this difference is maximally 0.5 and 0.3 dex 
for the volume- and flux-weighted metallicities, respectively. 
Because the spread in metallicity is much larger a.t z — 3 
than a.t z — 0.25, the flux-weighted metallicities trace the 
mass-weighted 84th percentile at 2 = 3 as opposed to signif- 
icantly higher metallicities at z — 0.25. 



6 EFFECT OF AGN FEEDBACK 

We repeated the surface brightness proflle analysis for sim- 
ulations including AGN feedback, because this feedback has 
been shown t o be important for the halo gas around massi ve 
galajcies (e.g. iMcCarthv et al.|[201(]| : iBertone et al. I l2010al lbl: 
Ivan de Voort fc Schavd \20m . Because the effect of AGN 
feedback is largest at low redshift and for soft X-ray emis- 
sion, we only show those results (Figure [7|, but describe the 
small differences for UV emission as well. 

The surface brightness proflles for soft X-ray lines are 
shown in Figure [7] as the dashed curves. The solid curves 



show the result for the simulations without AGN feedback 
and are identical to those shown in Figure [T] AGN feedback 
reduces the emission from 10^^"" Mq haloes by more than 
an order of magnitude for N vil, O vill, and Ne x in the 
halo centre. The reduction is smaller at larger radii, but 
still 0.5 dex at 0.3i?vir. For group-sized haloes, the emission 
in the central region is also lower by 1 dex for Ne x, 0.6 dex 
for N vil and O vill, 0.3 dex for C vi, whereas the O vil 
surface brightness in the centre remains the same. At larger 
radii, the reduction is small, about 0.2 — 0.3 dex at 0.5-Rvir. 
In the outer region of cluster haloes, there is no difference 
when we include AGN. However, in the most inner region, 
R < 0.1-Rvir, the surface brightness of all the soft X-ray lines 
increases by 0.1-0.7 dex when AGN feedback is included. 

The predictions for the maximum radius out to which 
the soft X-ray emission is detectable are not affected much 
by AGN feedback, except for the lowest halo masses. ORI- 
GIN and IXO/ATHENA can detect O viii out to 0.15i?vir, 
0.6i?vir, and 0.7i?vir for haloes with Mhaio = 10^^"", 
10"~", and 10""^^ Mq, respectively. C vi and O vil have 
surface brightnesses above 0.1 photon s~^ cm^^ sr^^ within 
0.15 — 0.3i?vir. N vil and Ne x drop below the detection 
threshold for Mhaio = 10^^"". 

The analysis for UV and Ha emission at z = 0.25 was 
repeated for a simulation including AGN feedback, but with 
lower resolution. We can therefore only compare the two 
highest mass bins. In the centres of haloes {R < 0.1-Rvir) 
with Mhaio = lO^^"^'^ M0 the surface brightness is reduced 
by 0.5-1 dex. The reduction is largest for Civ and Ovi 
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Figure 7. Surface brightness profiles as in Figure[T]for a simulation without AGN feedback (solid curves, identical to those in Figure[T]| 
and for a simulation with AGN feedback (dashed curves). 



emission. At larger radii, the reduction is only about 0.2 dex. 
Surface brightnesses are not affected by AGN feedback for 
Mhaio = 10^^"^" Mq at i? > 0.03i?vir and are only slightly 
reduced at smaller radii. 

AGN feedback affects the rest-frame UV surface bright- 
ness profiles at z = 3 the most for haloes with Afhaio = 
-|^gii-i2 where the reduction in the centre {R < 0.1i?vir) 
is 0.3 dex for Cm, Siiii, and Siiv and up to 0.8 dex for 
Civ and Ovi. At larger radii, the maximum reduction is 
0.2 dex. The implemented AGN feedback may be too strong 
for haloes of this mass for the relatively high resolution of 
the 25 h~^Mpc simulation used at z — 3, so the actual 
surface brightnesses could be closer to the ones predicted by 
the simulation without AGN feedback, as shown in Figure [21 
For haloes with A/haio ~ 10^""^^ M0, the surface brightness 
profiles are lower (by 0.2-0.5 dex) with AGN feedback, but 
for haloes with Mhaio ~ 10^^~^^ Mq, the surface brightness 
profiles of C lll. Si lll, and Si iv are actually higher with AGN 
feedback (by 0.2-0.4 dex). 

The mass-, volume-, and flux- weighted physical prop- 
erties we find in the simulations with AGN feedback (not 
shown) are similar to those in the simulations without AGN 
feedback. The differences are largest for X-ray emission. 
The weighted temperatures change by less than 0.1 dex. 
AGN feedback reduces the median flux-weighted overdensi- 
ties by 0.1-0.4 dex for most halo masses. Mass- and volume- 
weighted metallicities increase by 0.1-0.4 dex, but the flux- 
weighted metallicities increase less when AGN feedback is 
included. 



7 DISCUSSION AND CONCLUSIONS 

A large fraction of the gas in galaxy haloes, groups, 
and clusters has temperat ures T — 10*'^~^ K (e.g. 
Ivan de Voort fc Schave|[2012l ). At these temperatures, cool- 
ing is dominated by metal-line emissio n for metallicities 
Z > 0.1 Zq (e.g. IWiersma et all [2009al '). Observing these 
lines therefore constitutes an excellent route to studying the 
diffuse, warm-hot gas around galaxies, both at low and at 
high redshift. Additionally, metal-line emission provides us 
with clues as to how the circumgalactic gas was enriched 



and thus about the feedback process responsible for this en- 
richment. 

We used cosmological, hydrodynamical simulations 
from the OWLS project to quantify the surface brightness 
proflles of diffuse, circumgalactic gas, for both soft X-ray 
and rest-frame UV metal lines. This was done by stacking 
surface brightness maps centred on galaxies in the simula- 
tions. We only included emission from gas with densities 
wh < 0.1 cm~^, so our predictions should be considered to 
be lower limits to the actual surface brightness predicted by 
the models. The lines we considered are C vi (0.367 keV), 
N vii (0.500 keV), Ovii (0.561 keV), O viii (0.654 keV), and 
Ne X (1.021 keV) in the soft X-ray band and C in (977 A), 
C IV (1548 A), Si III (1207 A), Si iv (1294 A), and O vi 
(1032 A) in the rest-frame UV band. We discussed their 
detectability with current and future instruments and com- 
puted the flux-weighted physical properties of the gas. We 
also quantifled the effect of AGN feedback on our results. 

Proposed X-ray missions with detection limits of 
10^^ photon s~^ cm~'^ sr~^ can detect metal-line emission 
in galaxy haloes, groups, and clusters at 2: = 0.125. Ovill 
emission is the strongest and can be detected out to 80 per 
cent of the virial radius of groups (Mhaio = 10^^~^^ ^0) 
and clusters (Mhaio = lO""^'^ Mq) and out to 0.4_Rvir for 
Mhaio = 10^^"^^ Mq. Cvi, Nvii, Ovii, and Ne x can be 
detected out to smaller radii, 0.1 — 0.5-Rvir. 

Assuming a detection limit of 

10"^'' erg s~^ cm~^ arcsec"^ for UV metal lines at z = 0.25, 
C III can be detected out to 0.3i?vir for 10^^"^^ Mq haloes 
and out to 0.5-Rvir for 10^^~^* Mq haloes. It is the strongest 
UV metal line, but because it has a wavelength blueward 
of Lya, it may be absorbed by intervening hydrogen 
clouds. Civ, Ovi, Sim, and Siiv can be detected out to 
0.1 — 0.2 _Rvir for halo masses above 10^^ Mq (the lowest 
mass we considered at low redshift), and Sim can even be 
seen out to the same radius as C ill for 10^'^"^* Mq haloes. 

Assuming a detection threshold of 
10"^" erg s~^ cm~^ arcsec"'^, upcoming optical instru- 
ments should be able to detect the rest-frame UV metal 
hues Civ, Ovi, Sim, and Siiv out to O.lTivir at z = 3 for 
haloes more massive than 10^^ Mq. C ill is the strongest 
emission line and can be observed out to twice that distance. 
However, as we noted above, C ill emission may be strongly 
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absorbed by intervening neutral hydrogen, especially when 
it is emitted at high redshift. The other UV metal lines 
could therefore be stronger. The same instruments can also 
observe Hq at low redshift, which is a good probe of the 
cold (10'' K) halo gas. Ha may be detected out to 0.2, 0.3, 
and 0.5 7?vir for Mhaio = 10"-^^ 10^^"", and lO""" Mq, 
respectively. 

The flux-weighted mean temperatures of the gas probed 
by the soft X-ray lines increase with halo mass, but less 
steeply than the mass- and volume-weighted mean temper- 
atures. The emission is biased towards high density gas, 
especially in the inner halo and for cluster-sized haloes. 
The volume- and mass-weighted mean metallicities are sim- 
ilar, but always a factor of 3-10 lower than the flux- 
weighted mean metallicities. The differences between the 
flux-weighted properties of different metal lines are small 
for soft X-ray lines, although for clusters Ovi probes sig- 
nificantly lower temperatures and higher densities than the 
other metal lines that we consider. 

The flux-weighted temperatures of the gas from which 
the rest-frame UV emission originates is very close to the 
temperatures at which the corresponding emissivity curves 
peak, and they are insensitive to halo mass. The temperature 
of the gas probed is similar (T ~ 10^'^ — 10^ K) for all lines, 
except for O vi, whose emissivity peaks at T ~ 10^'^ K. 
At low redshift, all the metal lines are biased towards high 
densities. This behaviour is strongest for the silicon lines and 
weakest for O vi. At high redshift, the bias is smaller. In the 
outer parts of haloes with Mhaio < 10^^ Mq , C iv and O vi 
have flux-weighted mean densities below the mass-weighted 
mean density, but still above the volume-weighted average. 
The volume- and mass-weighted metallicities are always a 
factor of 3-30 lower than the flux-weighted values. 

AGN feedback introduces significant uncertainties in 
the predicted surface brightness profiles, especially at small 
radii. With AGN feedback, the emission decreases for low- 
mass haloes, but increases for some of the highest-mass 
haloes. In most cases, the radius out to which we expect 
to observe the metal-line emission is, however, not affected 
much. 

In this work we have not tried to mimic observational 
errors and selections effects and we have not accounted for 
some potentially important physical processes, such as ab- 
sorption by intervening clouds and emission from the multi- 
phase interstellar medium. Ifence, although our results can 
help build intuition about the level and interpretation of 
metal-line emission from diffuse gas in haloes, our predic- 
tions for the detectability of specific metal lines must be 
considered highly approximate. Future work could try to im- 
prove the realism of the predictions for specific instruments 
by creating virtual observations. 

Future work could also investigate the usefulness of 
more sophisticated statistical analyses. For example, stack- 
ing randomly oriented haloes, as we have done here, will not 
show filamentary emission, because the filaments will be ori- 
ented randomly. However, before stacking the galaxies, one 
can first rotate them towards their nearest neighbour of sim- 
ilar mass. This should enhance the emission of the gas in one 
direction. The signal at large radii is hard to detect, so the 
brightest emission lines should be chosen for this. We carried 
out some preliminary tests of this strategy using our simula- 
tions, but found that the enhancement is only about 0.2 dex 



(not shown). It may, however, be worthwhile to investigate 
other ways to rotate the galaxies before stacking. 

Besides stacking all images, it is possible to stack only 
pixels with a detection in the brightest emission line and 
then look for emission from other lines. Because the emis- 
sion is highest in the densest regions, the emission from 
different lines will be correlated. This will also lead to an 
enhancement of the signal. Pixels with a detection in Lya 
could, for example, be well-suited for detecting C iii, C iv. 
Si III, and Si iv as these metal-lines probe relatively cold 
(T ~ 10'' — 10^ K) gas. We intent to investigate such a 
strategy in future work. 
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APPENDIX A: RESOLUTION TESTS 

We tested the dependence of our results on pixel size, bin 
size of the surface brightness profile, and on the thickness of 
the region for which the emission is added. We found that 
the pixel size and bin size are unimportant for our results. If 
the thickness of the region is increased by up to an order of 
magnitude, the surface brightness outside the virial radius 
is increased by at most 0.3 dex. These surface brightnesses 
are below the unresolved X-ray background and therefore 
not observable with proposed instruments. 

Our results do not depend on the box size of the simu- 
lation, but there is a small dependence on resolution, which 
is shown in Figure lAll Blue curves show the original sur- 
face brightness profiles for the brightest line in Figures [T] 
[21 and O Red curves show the results for simulations with 
eight (two) times lower mass (spatial) resolution. Differences 
for the other X-ray and UV lines are very similar and thus 
not shown. The largest discrepancies between the low- and 
high-resolution simulation are found at small radii. 
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Figure Al. Mean (solid curves) and median (dotted curves) O VIII surface brightness at 2 = 0.125 (top row), C III surface brightness 
ai z = 0.25 (middle row), and C in surface brightness at 2 = 3 (bottom row) as a function of radius for the same haloes as in Figures[T] 
[2] and [S] respectively. Note that different rows show different halo mass bins. The blue curves show the profiles for high-resolution 
simulations, as also used throughout this paper. The red curves show the profiles for simulations with eight (two) times lower mass 
(spatial) resolution. 
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